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.4 bstract 

An Hxisv'tnini'rrir traveling' magnetic wave induces a meridional base flow in a cylindrical zone of 
an electrically conducting liquid. This remotely induced flow can be conveniently controlled, in 
magnitude and direction, and can have benefits for crystal growth applications. In particular, it 
can lie vised to significantly offset natural convection. Theoretical basics of this new technological 
method are presented. 


1 . Introduction 

Buoyancy convection, usually present 
during single crystal growth, is known to 
play a crucial role in defect formation 
and quality of the crystal. It originates 
from externally imposed thermal and/or 
solutal gradients necessary for 
crystallization, and the presence of 
gravity. Considerable research has been 
focused on controlling convection with 
the aim of improving macro- and 
microscopic homogeneity of the grown 
material. Utilizing a microgravity 
environment is the obvious solution to 
significantly reduce the gravity level and 
hence convection. Another method 
involves the application of static or 
dynamic magnetic fields. This method 
can be used only for electrically 
conducting melts, such as liquid metals 
and semiconductors, as well as ionic 
solutions or melts. While static magnetic 
fields can only play a passive role to 
damp already existing flows, time 
dependent magnetic fields can set up a 
steady state Lorentz body force, 
inducing flow in the system. The most 
recognized example is the rotating 
magnetic field (RMF). Over the past 
decade, RMF has been extensively 
studied as means for generating 
convective patterns beneficial to 
solidification processing [1,2]. It has been 
used for such growth techniques as 
vertical Bridgman( VB), float zone(FZ), 


or the traveling heater method (THM). 
The molten zone in these techniques 
employs a cylindrical geometry. When 
RMF is applied to such a zone, it 
induces a swirling, rotating, flow. This 
flow, in turn, builds up a secondary flow 
in a form of two rolls, one adjacent to 
the bottom of the molten cavity, and 
the other at the top. The secondary flow 
direction is always pointing into the 
melt at the central axis. The benefits of 
this type of flow pattern are: a) the 
azimuthal flow homogenizes and 
reinforces axisymmetry in the melt, b) it 
may suppress oscillatory buoyancy 
convection. Crystals grown in presence 
of RMF are reportedly of better quality 
[3-5]. We have recently proposed yet 
another type of magnetic field 
configuration that may find use in 
crystal growth practices[6,7]. It consists 
of an axisymmetric traveling magnetic 
wave that induces a steady state 
Lorentz force in the axial direction. 
When applied to a long open tube, it 
will pump the fluid. For closed cavities, 
a recirculating meridional flow is set up. 
These are essential features of the 
traveling magnetic field (TMF) method. 
Consider now a VB configuration with 
moderate radial thermal gradients 
(Grashof number Gv of the order of 10 5 
or less). Natural convection will then 
consist of a single axisymmetric roll. 
Application of TMF to this system w T ill 
significantly modify the flow pattern. 
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convection. This possibility may be 
at t rant i vo fo r in tor face shape control, 
reducing concavity, for example. A 
practical implementation of the 
technique consists of a stack of windings 
around the cylinder. A controlled phase 
shift of the applied electric power to 
each of these segments can result in a 
propagating magnetic wave. The power 
level, the frequency, the spacing between 
the sections and direction of propagation 
are the control parameters of this 
technique. Further feasible modifications 
include axially non-uniform TMF, TMF 
in combination with the static magnetic 
fields, TMF plus RMF, etc. 

Possible materials processing appli- 
cations of TMF include VB, FZ and 
THM crystal growth methods. TMF can 
be useful also in homogenizing melts in 
long cylindrical ampoules prior to 
solidification, an especially useful feature 
when processing in a microgravity 
environment. Finally, it can be used to 
control the sedimentation process taking 
place in ionic solutions. 


(VuisiiUT an infmiiely long cylinder of 
radius u, and electrical conductivity a. 
We further assume that the surrounding 
cylinder medium is an insulator and a 
11011 -magnet ic material. An external 
axisy m metric magnetic wave is applied 
to this cylinder along itvS axis. The 
frequency of this field is such that no 
propagating electromagnetic field is 
generated in the surrounding cylinder 
space. In this particular case, the 
magnetic field can be conveniently 
represented by its vector potential A, as 
it has only an azimuthal component. A 
single harmonic w T ave propagating in the 
z-direction along the axis of the cylinder 
of the frequency uj and the wave-vector 
a can be represented as: 

A = eA{r)e l ^- iU> (1) 

Here we follow the notation that the real 
parts of the complex variables 
correspond to physical quantities. The 
following well-known set of 
electromagnetic equations is employed in 
our model: 

B = V x A , VxH = j 


In this work we examine the possibility 
of counter-balancing natural convection 
by application of TMF for VB and THM 
configurations. In order to demonstrate 
damping of the natural convection 
effect, we consider here a case of VB 
with germanium in a silica ampoule. We 
numerically obtain the buoyancy force 
distribution and then consider a TMF 
induced force superimposed upon it. The 
combined force can be, for example, 
optimized is such a way that the flow in 
the vicinity of the growing interface is 
minimized. Also, a model of THM with 
applied TMF is presented and a 
buoyancy convective flow reduction as 
much as 50 times is demonstrated. 


VxE = - 
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where B, E and <f> are the magnetic 
inductance, electric field, and electric 
potential, respectively. It is convenient 
to select the Coulomb gauge for A: 
VA = 0 . The electric current density 
induced in the electrical conductor can 
be written in the following general form: 


j = <t -V0 


DA 

dt 


+ v x B 


(3) 


where v is the local velocity of the 
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and using the selected gauge for A we 
for the electric potential V 2 b — 0 . 
As both j and A have only the 
azimuthal components, then components 
of these vectors normal to the interface 
are automatically zero. Consequently, 
the derivative of o normal to the 
interface is also zero, and in such a case 
the solution for 0 is const , which can be 
set to zero. Substituting a general 
axisymmetric form (1) into Eq.(4) 
results in the following differential 
equation: 


d 2 Id 
dr 2 r dr 
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A = 0 


( 5 ) 


with the boundary condition of 
continuity for A together with its 
normal derivatives across the interface. 


It is convenient to introduce a non- 
dimensional scaling of the distance. The 
most natural scale in our model is the 
cylinder radius r 0 and it is adopted for 
the rest of the paper. The resultant 
differential equation (5) can be easily 
integrated and leads to the following 
solution: 
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This field induces a net Lorentz force in 
the stationary medium with the steady- 
state component in the 2 -direction: 


. (s) 

For the low frequency regime, K<<L 
we obtain {F r j = 0 and 


/t-i\ u/’era a 2 t 2 

(C) = — A J^ ar \ (9) 

Note that this low frequency limit 
formula is valid for a finite cylinder 
geometry as well. In this limit, the 
induced magnetic field is neglected. For 
an infinite geometry, as frequency 
increases, the r -component of the force 
becomes noticeable and squeezes the 
cylinder. For a finite cylinder geometry, 
the numerical solution of the induced 
Lorentz force indicates significant cap 
effects for values of K in the range of 3- 
10 [9]. In this work, we confine analysis 
to low frequencies, K< 1, which is often 
the case for practical purposes. 


5. Numerical work 


The TMF concept emerged from 
recognizing deficiencies in the RMF 
technique to generate an optimal flow 
for crystal growth in cylindrical 
ampoules. RMF can be viewed as a 
magnetic wave propagating in the 
azimuthal direction along a closed path. 
The obvious requirement here is that 
admissible values for the wave- vector 
are integers; the other geometries 
involve open paths. The waves can be 
planar, axisymmetric, or more complex. 
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Fig. 1 . The anemometer response to a 
downward traveling magnetic field for 
three different values of applied 
magnetic field. 

However, the axisymmetric case is 
particularly attractive as it can be easily 
obtained and has potential applications. 

In order to demonstrate the basic 
principles of this novel technique, we 
conducted an experiment with a 
cylindrical column of mercury [6], The 
radius of the column is 3.175 cm and its 
height is 12 cm. A stack of six coils was 
wound around the column with a 2.5 cm 
length of a single section. The frequency 
of the applied voltage was 50 Hz. The 
flow response to the TMF strength was 
studied by utilizing a hot film 
anemometer. Fig.l represents the 
anemometer signal evolution as the 
applied TMF power increases stepwise 
from zero. The rms value of the applied 
magnetic field was measured at the 
surface of the column. A clear transition 
from the steady to oscillatory flow 
occurs between 3.5 and 7.2 Gauss. 
Incidentally, this value is close to the 
unsteady flow transition for the RMF 
technique. Application of TMF power 
above 14 Gauss results in turbulent 
motion. A millitesla level of the 
magnetic field required to induce a 
significant flow is an advantage of this 
technique. The induced Lorentz force is 
in the range of the buoyancy force 



Fig. 2. The geometry and thermal fi<dd 
contours of the VB and THM cells used 
for numerical modeling. 


commonly encountered in solidification 
processing. It becomes then apparent 

that the TMF can potentially 
counterbalance the buoyancy force. A 
stagnant layer can be formed at the 
growing interface, or the flow direction 
can be reversed. Clearly, accurate 
balancing the two forces can be 

experimentally challenging; a flow 

transducer may be required for this 

purpose. Flow reversal can be 
accomplished with less difficulty. 
Experimental work in this direction has 
just begun, and our goal here is to 
numerically demonstrate the feasibility 
of damping convection by TMF. VB 
and THM models are selected for this 
mostly illustrative purpose. The analysis 
is further restricted to the steady-state 
axisymmetric case. Also, the low 
frequency approximation (9) is adopted 
for this calculation. The Lorentz force 
averaged over one oscillation for the 
non-zero motion of the medium contains 
terms linear in velocity which are not 
present in Eq.(9). With these terms 
included, the following equation for the 
vorticity, induced by a combination of 
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rhe heat rranspurt. equation needed tor 
buoy aiiey convection analysis. 

Additional, not explicitly displayed, 
equations are the incompressibility 
condition and the definition of the 
vorticity. In general, we have a multi- 
parametric problem described by a set of 
non-dimensional parameters ; Pr , Gr ? 
Tm, Ha, L, a, K. We confine our 


where the adopted here usual notation 
and scaling is: 
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where T top is the temperature at the top 
of the cell, T boti is the temperature at the 
bottom of the VB cell or at the center of 
the side wall for the THM cell, Pr is the 
Prandtl number, Tm is the magnetic 


analysis here to the case of P/ — 0.1. a=l. 
K< 1, Ha< 1, and the aspect ratio L= 1. 
Fig. 2 represents the geometry and 
thermal profiles for both the V B and 
THM models. The VB model describes 
molten germanium above a germanium 
crystal enclosed in a quartz ampoule. 
The temperature along the external wall 
of the ampoule is assumed to vary 
linearly and the top and bottom surfaces 
are isotherms. The ratio of thermal 
conductivities of the three materials, 
molten Ge, solid Ge, and silica is 
51:16.7:1.5. correspondingly. Due to this 
significant difference in thermal 
conductivities, the radial thermal 
gradient sets up at the growing crystal 
interface and induces natural 
convection. For Grashof numbers less 



r 

Fig. 3. A stream function contour plot of meridional flow induced in the VB cell by a 
combination of buoyancy and TMF force. The TMF applied force increases from the 
left to the right. At the 5 th and 6 th plots a stagnant fluid layer is formed at the bottom 
of the cell. 
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Fig. 4. Profiles of the r-componeut of 
flow along the z axis at r=0.6 for VB 
configuration and for 30 values of Tm. 


Fig. 5. Profiles of the z-cornponent of 
flow along the radius of the THM 
cell at z— 0 for a series of Tm values. 


than 10 6 . convection is steady and 
laminar[8]. The driving convection term 
in the vorticity equation (10a) is 
significantly non-uniform in the z - 
direction and it is confined to the 
growth interface area. In order to 
effectively counter-balance this term, the 
TMF term also has to be non-uniform. 
In this work a simplified case of the 
TMF, uniform in the z-direction. is 
considered. A commercial finite element 
with cubic basis numerical solver has 
been utilized for this work and adaptive 
gridding technique with approximately 
10000 nodes is implemented. 

When power of the TMF is applied in a 
direction opposite to the gravity force, it 
creates a small vortex at the top of the 
molten zone. At the same time, the flow 
in the buoyancy-induced vortex 
diminishes. As the TMF strength 
increases, the top vortex expands, until 
a stagnant layer is formed at the 
bottom. Further increase of the TMF 
field will set up a single vortex with the 
flow in a reverse direction to that of the 
initial natural convective flow. A two 
times reduction of the maximum flow 
velocity in the cell is obtained for some 
intermediate value of Tm. The stream 
function evolution of this process is 


depicted in Fig. 3 and the r -component 
of the flow as a function of z for r=0.6 
and for a set of Tm values is presented 
in Fig. 4. Clearly, a stagnant layer at the 
growth interface is observed. However, 
only for a narrow range of the magnetic 
Taylor numbers can its presence be 
detected. 

Our THM model consists of a central 
molten zone between the two solids 
encased into the cylindrical ampoule as 
depicted in Fig. 2. For the sake of 
simplicity, all three thermal 
conductivities are taken equal. At the 
side wall, a gaussian thermal profile is 
imposed and both the top and bottom 
surfaces are isotherms (see Fig. 2). The 
radial thermal gradient sets up a natural 
convection in the cell in a form of a 
single roll. As the vertical dimension of 
the molten zone is significantly lower 
than in the VB case, counter-balance of 
the buoyancy force by the uniformity in 
z-direction TMF can be more successful. 
Indeed, a significant reduction of the 
global flow is obtained. Fig. 5 represents 
a set of profiles of the z-component of 
flow velocity along the horizontal central 
line (z=0, r) for an equidistant sequence 
of Tm values. A damping of original 
flow as much as 50 times is obtained. 
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Summary 

In this work wi‘ have numerically 
demonstrated the feasibility of 
significantly offsetting gravity induced 
flows by the application of axisymmetric 
traveling magnetic fields. For melts with 
small aspect ratios, the global flow can 
be dramatically reduced, as demon- 
strated here, by 50- fold with THM 
application. For long cylinders 
convective flows induced by localized 
bony ancy forces arc not effectively 
reduced by application of uniform TMF. 
However, a stagnant layer of fluid 
adjacent to the crystal growth interface 
can be obtained. A better approach is to 
consider a non-uniform localized TMF. 
In fact, only four coils are sufficient to 
set up a traveling magnetic field. Finally 
we want to note that in many instances, 
a steady state single axisymmetric 
vortex flow is beneficial during 
solidification processing. In such cases 
TMF can be a solution. 
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